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Abstract 

We have measured the inchisivc semileptonic branching fractions of D^, D^, and mesons. 
For these measurements, we have used the fuh CLEO-c open-charm data samples, 818 pb^^ at 
Ecu = 3.774 GeV, giving D^D'^ and D+D- events, and 602 pb"^ at Eqm = 4.170 GeV, giving 
D*'^Df events. We obtain ^ Xe+u^) = (6.46 ± 0.09 ± 0.11)%, B{D+ Xe+u^) = (16.13 ± 

0.10 ± 0.29)%, and i3(L>+ Xe+i/g) = (6.52 ± 0.39 ± 0.15)%, where the first uncertainties are 
statistical and the second are systematic. From these and lifetimes obtained elsewhere, we obtain 
the ratios of semileptonic decay widths r(L>+ Xe+z/e)/r(L>o Xe+Ue) = 0.985 ± 0.015 ± 0.024 
and r{D+ Xe+Ue)/r{D'^ Xe+Ve) = 0.828 ± 0.051 ± 0.025. The ratio of D+ and £)° is 
consistent with the isospin symmetry prediction of unity, and the ratio of Df and differs from 
unity, as expected. 
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I. INTRODUCTION 



As part of the CLEO-c analyses of exclusive [l|-[7|] and inclusive semileptonic decays |8|], 
this article presents measurements of D^, D~^, and inclusive semileptonic branching 
fractions using the complete CLEO-c data sets. Using these results and known lifetimes, we 
also report the ratios of the widths T{D^ Xe"*"z/e)/r(D° — )■ Xe'^i'e) (which is expected 
to be unity due to isospin symmetry) and T{D^ — Xe^i'e)/T{D^ — )■ Xe^u^) (which is not 
expected to be unity jo], [lo[, though with poor theoretical precision). These measurements 
are important in their own right, and, due to similarities between the D and B sectors, 
will also improve understanding of B semileptonic decays. In particular, knowledge of the 
previously unmeasured ratio T{D^ — )■ Xe^i>e)/r{D^ — )■ Xe'^i'e) enables a more reliable 
prediction of the difference of the inclusive decay rates between B^ and B~^ mesons in 
b — i- uC^vi decays, thereby reducing theoretical uncertainty (of in determination of weak 
mixing parameter Vuh- 

Two sets of open-charm data samples are used to study the semileptonic decays of charm 
and charmed-strange mesons. In e^e~ collisions provided by the Cornell Electron Storage 
Ring (CESR), the CLEO-c detector has collected integrated luminosities of 818 pb~^ at 
the center-of-mass energy Ecu = 3.774 GeV near the peak of the '?/'(3770) resonance which 
decays to DD pairs, and 602 pb~^ at Ecu = 4.170 GeV near the peak production of Dg'^Df 
pairs. The former data set contains 3.0 x 10® D^D^ and 2.4 x 10® D^D~ pairs, and is used to 
study and semileptonic decays. The latter data set contains 0.6 x 10® D*g^DJ pairs, 
and is used to study semileptonic decays. We have previously reported [1] measurements 
of inclusive semileptonic decay branching fractions of and mesons with a subsample 
of the former data set. 

The remainder of this article is organized as follows. The CLEO-c detector is described 
in Sec. [Tll Event reconstruction and selection criteria are described in Sec. Illli The analysis 
procedure to extract semileptonic decay rates is covered in Sec. [IVl Results for inclusive 
spectra are presented in Sec. |Vl Systematic uncertainty in our measurements is evaluated 
in Sec. |Vll Finally, in Sec. IVIIl a summary of our results is provided. 



II. THE CLEO-c DETECTOR 



The CLEO-c detector [TMi is a general-purpose solenoidal detector equipped with four 
concentric components: a six-layer vertex drift chamber, a 47-layer main drift chamber, a 
ring-imaging Cherenkov (RICH) detector, and a cesium iodide electromagnetic calorimeter. 
The detector provides acceptance of 93% of the full 47r solid angle for both charged particles 
and photons. The main drift chamber provides specific- ionization (dE/dx) measurements 
that discriminate between charged pions and kaons. The RICH detector covers approxi- 
mately 80% of An and provides additional separation of pions and kaons at high momentum 
(> 700 MeV). Electron identification is based on a likelihood variable that combines the in- 
formation from the RICH detector, dE/dx, and the ratio of electromagnetic shower energy 
to track momentum (E/p). A GEANT-based 15|] Monte Carlo (MC) simulation is used to 
study efficiencies of signal and background events. Physics events are generated by evt- 
GEN f]^, tuned with improved knowledge of charm decays [l7-20], and final-state radiation 
(FSR) is modeled by photos [21|]. 
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III. EVENT SELECTION 



Charm or charmed-strange mesons are always produced in pairs in our open-charm data 
samples. Since the data are taken just above threshold, the mesons are produced in a very 
clean environment with no additional particles except, in the case of the DgD*, a photon or 
a neutral pion from the D* decay. The analysis proceeds by first defining a single tag (ST) 
sample, in which one of the D (or Dg) mesons in a DD (or DgD*) event is reconstructed 
in a chosen hadronic decay mode, and a further double tag (DT) subsample in which an 
additional recoiling electron (or positron) is required as a signature of the signal semileptonic 
decay. Absolute semileptonic branching fractions for charm or charmed-strange mesons can 
then be obtained from the fraction of the ST sample that is DT, without requiring any 
knowledge of the integrated luminosity or how many mesons are produced. 



A. Tag Selection 

To minimize the combinatorial backgrounds and systematic uncertainties, three very clean 
tag modes composed of only charged particles are used: — > K^ti^, — t- K^ti^ti^, and 
D~ — )■ 07r~. Here, the notation Dj — )■ 07r~ is a shorthand label for Dj — t- K^K^tt^ events 
within a 10 MeV mass window of the meson peak in K~K~^ invariant mass. The inclusion 
of charge conjugate modes is implied throughout this article unless otherwise stated. 

We identify a ST in the '?/'(3770) data sample using the energy difference AE = -Ed— -Ebeam 
and the beam-constrained mass difference AMbc = [-E^beam ~ Pz)]^''^ — ^d, where Ed is the 
energy of the tag, Eheam is the beam energy, p_D is the three momentum of the tag, and mo is 



the nominal mass 17| of the neutral or charged charm meson. We require the — )■ K^ir' 
and — )■ K^n^ii^ tags to have AMbc within a 4 MeV mass window around the nominal 
D mass. 

For data collected at the center-of-mass energy of 4170 MeV, we identify a ST by using 
the invariant mass of the tag M{Ds) and recoil mass against the tag Mi.ecoii(-Ds). The recoil 
mass is defined as M^ecoii{Ds) = [(-^ee - E^.Y - (Pee - PdJ^]^^^, where {Eee, Pee) is the net 
four-momentum of the e+e~ beam taking the finite beam crossing angle into account, and 
{Eds, PdJ is the four- momentum of the tag, with Ed^ computed from p^^ and the nominal 
mass llTi of the Ds meson. We require the recoil mass to be within 55 MeV of the D* 
mass [17|. This loose window allows both primary and secondary (from D*^ or 
D*^ —7- D^TT^) Ds tags to be selected. We veto tag candidates with track momenta below 
100 MeV to reduce the background from DD* decays (through D* — )■ ttD). 

The AE and AM distributions obtained from data are shown in Fig. [T] To estimate the 
backgrounds from the wrong tag combinations, we use the sidebands of the AE distribution 
or the tag mass difference AM = M{Ds) — distribution, where m^^ is the nominal 



mass 



17j of the Dg meson. We define the signal and sideband regions in Table HI We fit 
the distributions to a sum of a double- Gaussian function (for signal) and a second order 
Chebyshev polynomial function (for background) to determine the tag sideband scaling 
factor Stag, which is the ratio of areas in the signal and sideband regions described by the 
background polynomial function. Obtained ST yields and tag sideband scaling factors are 
listed in Table [TTl 
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TABLE I: Signal and sideband regions of AE and AM for each tag mode. 



Tag mode 



Signal (MeV) 



Sideband (MeV) 



-30 <AE< +30 
-25 < < +25 
-20 < AM < +20 



-80 <AE< -50 
+50 <AE< +80 
-65 <AE< -40 
+m <AE< +65 
-55 < AM < -35 
+35 < AM < +55 



1851209-010 




CD 

5000 



40000 



30000 



20000 



- 10000 - 



1 


1 


: D" ^ 71" 71" 

- 




J 





2000 



1500 - 



1000 - 



- 500 - 




-0.1 -0.05 0.05 0.1 -0.1 -0.05 0.05 0.1 

AE (GeV) AE (GeV) 



-0.05 0.05 

AM (GeV) 



FIG. 1: Tag AE and AM distributions in data (histograms) with fits (solid curves) and background 
contributions (dashed lines). 



B. Signal Selection 



We form DT candidates from ST candidates by adding a recoiling charged track that is 
consistent with coming from the nominal interaction point. Specifically, the recoiling track's 
point of closest approach to the origin must be within 5 cm of the interaction point along 
the beam line and within 5 mm of the interaction point in the plane transverse to the beam 
line. We require the momentum of the track to be p > 200 MeV and the angle with respect 
to the beam to be | cos6'| < 0.80 so that all charged-particle identification (PID) information 
{dE/dx, RICH, and E/p) is available. The signal track in the DT candidate is also required 
to be identified as an electron, a charged pion, or a charged kaon, for further analysis. This 
is discussed in the next section. 



TABLE II: ST yields and statistical uncertainties in data, where ngrp is the yield in the tag signal 
region, n^rp is the yield in the tag sideband region, stag is the tag sideband scaling factor obtained 
from a fit to tag AE (or AM) distribution, and n^T is the scaled sideband subtracted ST yield. 



Tag 


mode 


"-ST 


"■ST 


Stag 


nsT 


^0 




144260 


2258 


1.067 


141851 + 383 


D- 




231429 


7748 


1.104 


222872 ± 490 


D7 


—?■ (j)Tr~ 


10453 


807 


0.979 


9663 + 106 
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IV. ANALYSIS 



The D (or Ds) semileptonic inclusive spectrum (or differential decay rate) can be ex- 
pressed as 

c^^SL _ 1 A^e _ 1 AripT/esL 
dp riD Ap nsT Ap ' 

where ud is the number of D mesons produced, is the number of produced primary 
electrons in bins of momentum p, nsT is the number of ST, An^T is the electron candi- 
date yield in bins of momentum, and esL is the (momentum-dependent) electron detection 
efficiency. The D semileptonic branching fraction can be obtained by integrating the dif- 
ferential spectrum and correcting for the 200 MeV momentum cutoff by extrapolating the 
spectrum below the cutoff. If we had a perfect MC modeling of the semileptonic decays, 
a simple momentum bin- by-bin correction factor could be used for esL- Instead, we use a 
more general unfolding [22] approach to minimize MC model dependence. 

The observed laboratory momentum spectrum ?/(6, ^track) of a particle identified as type 
h (= e, TT, or K) in bins of measured track momentum bin Ztrack can be modeled as a folded 
distribution. It is related to the true laboratory momentum n{a,j) via detector-response 
matrices that account for resolution and efficiency: 

y{b, ^track) = ^ ApiB{b\a, itrack) ^ Arack (^track | J )n(a, j) , (2) 
a j 

where a (= e, fi, tt, or K) is the true particle species index, n{a,j) is the true laboratory 
momentum spectrum in bins of true laboratory momentum bin index j of a particle type 
o, Arack(^track|a, j) is the tracking efficiency matrix, which describes the probability of a 
particle of type a with momentum in bin j to be reconstructed in track momentum bin 
^track, and ApiD(^|a, itrack) is the PID efficiency matrix, which describes the probability of a 
particle of type a with measured momentum in bin Ztrack to be identified as PID type b. We 
unfold j22| Eq. ([2]) to obtain the true momentum spectrum 



n(a 



e,JJ 



^tra 



track l^track 1 0- — 6, 



^track 



E 



Aplj^{b\a, Ztrack)l/(&, ^track) 



(3) 



where the A~^^s are the unfolded inverses of each efficiency matrix. Because we are interested 
in the primary electron laboratory momentum spectrum (to obtain the branching fraction) 
we use the electron solution after PID unfolding (a = e). 

In addition to finite resolution and efficiency, modeled by detector-response matrices, we 
have to consider possible backgrounds in our observed spectrum. We remove combinatorial 
wrong-tag background contribution by AE (or AM) sideband subtraction. Charge sym- 
metric nonprimary true electron backgrounds (from 7 conversion and tt^ Dalitz decay) are 
subtracted by using the wrong-sign (WS, opposite to the expected primary electron charge) 
electron sample. In the following subsections, we break the analysis described above into 
discrete steps. 
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A. PID Yield 



From a set of signal candidate tracks, we measure the PID yield y{b, i) in bins of PID type 
6, track momentum bin Ztrack? (or AM) signal and sideband regions isB, and right-sign 
(RS) or wrong-sign (WS) bin depending on the charge of the track and the flavor of 
the tag, where i is a collective index for (itrack, ^sb, "^rw)- The charge of the daughter kaon 
defines the flavor of the — )■ K^t[~ tag, and the charge of the tag defines D~ — )■ K^i[~t:~ 
and D~ — )■ 07r~ tags. The RS track is defined to be the track with the same charge as the 
tagged -D° daughter kaon or to be the opposite charge of the charged tags, and the WS track 
is defined the other way around. 

B. PID Unfolding 

We correct for PID efficiency and mis-PID crossfeed backgrounds using 

y{a, i) = Apl^{b\a, i)y{b, i), (4) 

where i is a collective index for (itrack, ^sb, "^rw)- The PID matrix ApiD{b\a) used in the unfold- 
ing is shown in Fig. |2l PID matrix elements associated with the charged pion are obtained 
from Kg — i- TT+TT^ events, the charged kaon elements are obtained from — )■ K^ti^ti^ 
events, and the electron elements are obtained from radiative Bhabha events (e^e~7) em- 
bedded in hadronic events. Here we treat muons as pions because muons in the momentum 
range in which we are interested behave almost the same as charged pions in the CLEO-c 
detector. The effect of this approximation is negligible on our branching fraction measure- 
ment because the probability of pions (and muons) to be misidentified as electrons is very 
small, as shown in Fig. [2J After solving the PID problem, we take the electron solution 
(a = e) for further analysis. 

C. Tag Sideband Subtraction 

To remove the wrong-tag combinatorial background, we perform AE (or AM) sideband 
subtractions after PID unfolding. After this process, we deal with real electrons from D (or 
Dg) meson decay. 

D. Wrong-Sign Electron Subtraction 

Charge symmetric secondary electrons are removed by subtracting the WS (secondary) 
electron yield from the RS (primary plus secondary) electron yield. After this process, we 
end up with primary electrons from D (or Ds) meson decay. 

E. Tracking Efficiency, ^track 

We obtain the tracking efficiency matrix Arack(^track| j) from MC simulation. This includes 
track finding efficiency and resolution effects. 
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FIG. 2: (color online). The components of the PID efficiency matrix ApiD{b\a) obtained from 
data. The matrix describes the probability of a particle of type o to be identified as a PID type b. 
We measured the PID matrix in momentum intervals of 50 MeV (some bins are wider due to low 
statistics) above the PID momentum cutoff 200 MeV. The cases with a ^ b, conventionally called 
the fake rate or mis-PID probability, are shown in points with statistical uncertainties. The cases 
with a = b, conventionally called the efficiency, are shown as solid lines. The discontinuities at 
momentum 700 MeV in fake rates and efficiencies are due to the fact that the RICH information 
is used for pion and kaon identifications only above 700 MeV. 



F. Tag Bias Correction 

The signal semileptonic efficiency esL requires a possible tag bias correction which would 
be introduced if the ST efficiency in the signal DT events is different from that when the 
other recoiling system is a generic D-meson (or D^-meson) decay. The effect of the tag bias 
can be express in terms of a ST efficiency ratio 



esL 



^bta 



(5) 



eST ^ST ^ST CsT <^ST 

where eoT is the DT efficiency, esT is the ST efficiency against generic decays in the recoiling 
system, egrp is the ST efficiency when the recoiling system is the signal semileptonic decays, 
ee is the signal electron detection efficiency given the tag in the other side is found, and 
6tag is a measure of tag bias in the efficiency Thus, fetag = Cst/^st and ee = eDx/e'sT- We 
expect this effect to be small due to chosen clean tag modes and low event multiplicity. 
We estimate tag biases in MC simulation: 6tag(^° ^ e+X) = 0.9965 ± 0.0017, 6tag(^+ ^ 

— )■ e~^X) = 1.0069 ± 0.0021, where uncertainties are 



1.0017 ±0.0021, and fetagP^ 



e+X) 

due to MC statistics. 



G. Doubly Cabibbo-suppressed Decay Correction 



Because of the doubly Cabibbo-suppressed decay (DCSD) and quantum correlation [23|, 
24i\ in coherent D^D^ production at the ?/'(3770) resonance energy, we need a correction 
for the observed semileptonic branching fraction using the — )■ K~7r~^ tag mode. The 
observed branching fraction i3obs requires a correction [2^ 

1 + i?ws 



24 



(6) 
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TABLE III: Summary of DT yields, statistical uncertainties, and correction procedure explained 
in Sec. llVi PID yields (Sec. IIV Ap for electron candidates (6 = e) are shown in the first group for 
tag signal region (S), tag sideband region (B), right-sign (R), and wrong-sign (W) bins, where 
the yields in the sideband region are scaled by the tag sideband scaling factor (Table [Tl]) for each 
tag mode. PID unfolded (Sec. IIV Bp electron yields (a = e) are shown in the second group. Tag 
sideband subtracted (Sec. lIV C]l electron yields are shown in the third group, followed by the wrong- 
sign subtracted yield (Sec. IIV Pp . tracking efficiency-corrected yield (Sec. IIV Ep . and remaining tag 
bias (Sec. HVFjl or DCSD (Sec. HVG]) corrected yield. 



Z)0 



D+ 



D- 



PID yield, electron candidates 

y{h = e,S, R) 

y{b = e,B,R) 

y{b = e,S,W) 

yib = e,B,W) 
PID unfolded yield, electrons 

y{a = e,S, R) 

y{a = e,B,R) 

y{a = e,S,W) 

y{a = e,B,W) 
Tag sideband subtracted electrons 

y{a = e,R) 

y{a = e,W) 
Wrong-sign subtracted electrons 
Tracking efficiency-corrected electrons 
Tag bias (and DCSD) corrected electrons 



6618.0 ±81.4 
41.6 ±6.7 

653.0 ± 25.6 

19.2 ±4.5 

7292.4 ± 90.7 
47.1 ± 7.7 
682.4 ±31.4 

21.3 ±5.3 

7245.3 ±91.0 

661.1 ±31.9 
6584.2 ± 96.4 

8361.0 ± 123.0 
8450.8 ± 124.3 



24834.0 ± 157.6 
332.4 ± 19.2 
711.0 ±26.7 

55.2 ±7.8 

27304.5 ± 174.8 
370.4 ±21.7 
812.8 ±33.8 
65.2 ±9.8 

26934.1 ± 176.2 
747.6 ± 35.2 

26186.5 ± 179.6 
33182.0 ± 228.2 

33125.6 ± 227.9 



553.0 ± 23.5 
24.5 ±4.9 
50.0 ±7.1 
9.8 ±3.1 

608.9 ± 26.4 
27.7 ±5.6 
56.7 ±8.6 
11.7 ±3.4 

581.2 ±27.0 
44.9 ±9.2 

536.3 ± 28.5 
681.3 ±36.4 
676.6 ± 36.2 



is the ratio of the DCSD rate to the Cabibbo-favored 
—7- K^TT^) is the ratio of the time- integrated 



Here = \ {K+n-\D^) / {K+n 
decay rate, and Rws = r(£'° ^ K+7r-)/T{D^ 
DCSD rate to the Cabibbo-favored decay rate. Using the world average 
we need a correction factor (1 ± i?ws)/(l — r'^) = [1 + (3.80 ± 0.05) x 
0.09) X 10-3] = 1.0072 ±0.0001. 
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values of these 
3]/[l - (3.35 ± 



V. RESULTS 

The final electron candidate yields are summarized in Table IIIII and efficiency-corrected 
laboratory momentum spectra are shown in Fig. [31 Also shown in Fig. 12] are the spectrum 
extrapolations below the PID momentum cutoff (200 MeV). The curves shown are obtained 
with a fit using the sum of measurements of exclusive channels together with form-factor 
models and adding higher-resonance and nonresonant channels to match the sum of the 
exclusive channels with our measured branching fraction. Further details of the extrapolation 
procedure are available in the Appendix. From the fit results, we obtain fractions below the 
momentum cutoff of 7.8% for 8.0% for D+, and 7.0% for Df. 

At this point, we also consider the secondary electrons from leptonic decays of — > r+z/^ 
and — t- r+z/^ as they produce electrons through r+ — e~*"i^e'^r decay. This source 
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FIG. 3: Inclusive laboratory frame electron spectra obtained from data, shown as points with 
statistical uncertainties. The vertical dashed lines indicate the PID momentum cutoff at 200 MeV. 
Extrapolated spectra are shown as solid curves. The dashed curve in the spectrum plot is the 
expected contribution from e^VeV-r from leptonic D'^ — t- t^Vt decay. 



TABLE IV: Summary of semileptonic branching fractions. Here ;Btrunc is the partial branching 
fraction above 200 MeV, B{e'^X) is the extrapolated full branching fraction, and B{Xe^Ue) is the 
semileptonic branching fraction after r — )• e correction (for D'^ and Df). First uncertainties are 



statistical and the second are systematic due to uncertainties in B(D^ — >• T'^Vfj,) [27|, B{Df 



T+Ur) I25|,l26|], and B{t 



Tag 


mode 


5trunc(e+X) (%) 


B{e+X) (%) 


B{Xe+Ve) (%) 


^0 




5.958 ± 0.084 


6.460 ± 0.091 


6.460 ± 0.091 






14.863 ± 0.092 


16.147 ±0.100 


16.129 ±0.100 ±0.000 




— 7- (pTi^ 


7.002 ± 0.361 


7.525 ±0.387 


6.522 ± 0.387 ± 0.079 



of secondary electrons is expected to be large in so we have included the expected 
spectrum component in the extrapolation. The expected branching fractions of these 
secondary electrons from the leptonic decays oi and are subtracted from the 
fully inclusive branching fraction results to obtain inclusive semileptonic decay branch- 



ing fractions. The branching fraction for t^Ut decay is taken from Refs. [25|, |26 

B{D^ —7- t^Ut-) = (5.62 ± 0.41 ± 0.16)%. The size of the expected secondary electron 
contribution from the unobserved leptonic decay — > t^Uj. is based on the known branch- 



ing fraction of — /i^z^^^ decay [27[ scaled by the standard model decay rate ratio [17 
r{D^ — 7- T^Ur)/T{D^ — 7- fJ'^i'fj,) = 2.67. We take the uncertainty in the r — )■ e correction as 
a part of our systematic uncertainty. Branching fraction results are summarized in Table IIVI 
with all above-mentioned efficiency and cutoff corrections. 

The laboratory frame electron momentum spectra shown in Fig. [3] are given in tabular 
form in Table |Vl 
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TABLE V: Inclusive semileptonic electron partial branching fractions of D , D^, and Df in the 
laboratory frame. For and D^, we have subtracted expected contributions from leptonic decays 
T^Vr (followed by e^z/e^v)- Systematic uncertainties in total branching fractions are added 

to the statistical uncertainties. In comparing theoretical predictions with these measurements, one 
must smear the theoretical predictions by boosting from the D (or Dg) rest frame to the laboratory 
frame. For Dg, 51% of the electrons are from secondary Dg from D*, and 49% are from primary 



Dg. 


p (GeV) 


AB{D^ Xe+Ue) (%) 


A T -\ / I "A _l_ X r _|_ \ / rw \ 

AB{D+ Xe+Ve) (%) 


A T— 1 / r ~\ _l_ T r _l_ \ / rt-/ \ 

AB{Df ^- Xe+Ue) (%) 


0.200-0.250 


0.347 ±0.036 


0.912 ±0.040 


0.491 ±0.152 


0.250-0.300 


0.426 ± 0.030 


1.133 ±0.038 


0.470 ±0.124 


0.300-0.350 


0.576 ± 0.031 


1.379 ± 0.041 


0.554 ± 0.126 


0.350-0.400 


0.629 ± 0.030 


1.462 ± 0.043 


0.515 ± 0.120 


0.400-0.450 


0.640 ± 0.031 


1.675 ± 0.047 


0.578 ±0.112 


0.450-0.500 


0.640 ± 0.031 


1.661 ± 0.046 


0.562 ±0.123 


0.500-0.550 


0.596 ± 0.029 


1.546 ± 0.044 


0.794 ±0.127 


0.550-0.600 


0.575 ± 0.029 


1.415 ±0.041 


0.611 ±0.115 


0.600-0.650 


0.492 ± 0.026 


1.243 ± 0.038 


0.471 ±0.104 


0.650-0.700 


0.374 ± 0.023 


0.946 ± 0.032 


0.314 ±0.087 


0.700-0.750 


0.269 ± 0.019 


0.674 ± 0.026 


0.246 ± 0.079 


0.750-0.800 


0.230 ±0.017 


0.429 ±0.019 


0.089 ± 0.060 


0.800-0.850 


0.089 ±0.011 


0.240 ± 0.014 


0.115 ±0.060 


0.850-0.900 


0.053 ±0.008 


0.103 ±0.009 


0.037 ± 0.046 


0.900-0.950 


0.021 ±0.005 


0.022 ± 0.004 


0.074 ± 0.051 


0.950-1.000 


0.002 ± 0.002 


0.004 ± 0.002 


0.096 ± 0.045 


1.000-1.050 






0.015 ± 0.022 



VI. SYSTEMATIC UNCERTAINTIES 

Possible sources of systematic uncertainties and their effects on the branching fraction 
measurements are summarized in Table IVI[ 

The ST yields are obtained from a tag {AE or AM) sideband subtraction method. 
Because of the chosen clean tag modes, there is very little combinatorial background under 
the signal peak, as shown in Fig. [H Systematic uncertainties in the numbers of tags are 
studied by using alternative signal and background functions, and comparing the known 
input number of ST in a MC simulation test to the output with the same procedure. By 
adding all of the resulting variations in quadrature, we obtain 0.5% (in Z)°), 0.7% (in -D^), 
and 0.9% (in Df) uncertainties in the estimation of the number of ST. 

The systematic uncertainty of 0.3% in tracking efficiency was estimated [l8| in a detailed 
MC and data efficiency comparison using iIj{3770) — )■ DD events with the cases when both 
D and D mesons can be fully reconstructed. 

Uncertainties in FSR and bremsstrahlung effects on D semileptonic decay branching 
fraction measurements were studied in our previous measurement |8|] and in high statistics 
exclusive D semileptonic decay modes [sj. They are found to be well simulated in our MC 
program. We have assessed the uncertainty in FSR by redoing the analysis using alternative 
signal efficiency and input spectra with FSR turned off in the MC simulation. Including the 
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TABLE VI: Summary of sources of systematic uncertainty and their effects on the semileptonic 
branching fraction measurements. 



Source 


(%) 


D+ (%) 


Dt (%) 


Number of tags 


0.5 


0.7 


0.9 


Tracking 


0.3 


0.3 


0.3 


PID 


0.8 


0.5 


0.6 


FSR 


0.5 


0.5 


0.5 


Tag bias 


0.2 


0.2 


0.3 


DCSD 


0.0 






T ^ e 




0.0 


1.2 


Extrapolation 


1.3 


1.4 


1.5 


Total 


1.7 


1.8 


2.3 



uncertainty in bremsstrahlung simulation [5|, we assign 0.5% uncertainty due to FSR and 
bremsstrahlung effects on our branching fraction measurements. 

Uncertainties in electron identification for semileptonic decays are assessed by comparing 
the efficiency measured using a radiative Bhabha sample embedded in hadronic events to 
those in various MC simulated event samples. We assign systematic uncertainties due to 
electron identification as 0.7% for — )■ Xe'^i'e, 0.5% for — )■ Xe'^i'e, and 0.6% for 
— )• Xe'^i'e. For other PID efficiencies, we have varied their values within measured 
uncertainties and observe the effect on our measured branching fraction. By adding all 
electron identification and other PID uncertainties in quadrature we assign uncertainties in 
PID as 0.8% for ^ Xe+u^, 0.5% for D+ Xe+u^, and 0.6% for Df Xe+u^. 

Tag bias corrections are estimated from MC simulation. We take the uncertainty in the 
MC statistics and a quarter of the size of the tag bias as the uncertainty in the correction. 

For the and inclusive electron spectra, we subtract the contribution from 

— 7- e'^Uei^T as estimated in Table IVlnl (and Table HXl) to obtain the inclusive semileptonic 
branching fraction. We take the uncertainty from B{D^ — )■ t+z/t-) [and B{D~^ — )■ t~^i/t-)] 
listed in the table for the uncertainty on the r — )■ e contribution correction. The uncertainty 
in — !■ Xe'^i'e is negligible, and we assign an uncertainty of 1.2% in — )■ Xe~^i>e. 

To estimate systematic uncertainties in the extrapolation procedure, we fix all parameters 
to the reference values listed in Table IVIIt Table IVIIIt and Table IIXI Then we vary each 
semileptonic decay component one-by-one within the allowed range of uncertainties, listed 
in the table, and reevaluate the fraction below the momentum cutoff and the effect on the 
resulting branching fraction. For the unobserved decay components, we vary 100% of the size 
of the predicted branching fraction to assess the uncertainty. We also use alternative form- 
factor models, by changing models component-by-component from the reference models in 
the tables, when we perform an extrapolation fit as described in the Appendix, to assess the 
additional uncertainty in the extrapolation. By adding all effects in quadrature, we assign 
1.3% for Xe+i/e, 1.4% for Xe+z/g, and 1.5% for Df Xe+z/g as uncertainties 

in the extrapolation procedure. 
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VII. SUMMARY 



Using the full sample of open-charm data collected by the CLEO-c detector, we obtain 
the charm and charmed-strange meson inclusive semileptonic branching fractions: 

^ Xe+z/e) = (6.46 ±0.09 ±0.11)%, 
^ Xe+z/e) = (16.13 ±0.10 ±0.29)%, 

and 

S(D+^Xe+z/e) = (6.52 ±0.39 ±0.15)%, 



where the first uncertainties are statistical and the second are systematic. Using known [17 
lifetimes Too = (410.1±1.5) x lO'^^ s, td+ = (1040±7) x lO'^^ s, and = (500±7) x lO'^^ 
s, we obtain the ratios of semileptonic decay widths 



r(DO ^ Xe+z/e 

and 



0.985 ±0.015 ±0.024 



0.828 ±0.051 ±0.025. 



In these ratios, we assume the PID and tracking uncertainties are fully correlated and all 
others are uncorrelated. The former ratio shows that charged and neutral charm meson 
semileptonic decay widths are consistent with isospin symmetry, as expected, because the 
two mesons differ only in the isospin of the light quark. On the other hand, the latter ratio 
shows that there is an indication of difference between charm and charmed-strange meson 
semileptonic decay widths. 
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Appendix: Spectrum Extrapolation 

Charm and charmed-strange exclusive semileptonic decay components used to perform 
spectrum extrapolation fits are summarized in Tables rvTT| I Villi and llXI Efficiency-corrected 
data points are fit to a sum of exclusive semileptonic decay components to estimate the 
unmeasured portion of the spectrum below the momentum cutoff at 200 MeV due to the 
electron identification. Normalization of each component is allowed to float within the 
uncertainty shown in the tables. 

Higher-resonance and nonresonant decay components are used to make the sum of exclu- 
sive branching fractions match the inclusive branching fraction in and extrapolations. 
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TABLE VII: Summary of semileptonic decays used to perform the spectrum extrapolation. 
Assumed branching fractions are shown in the second column; normalization of each component 
is allowed to float within the given uncertainty. Form-factor models used to describe the shape of 
each spectrum are shown in the third column: single-pole (SPOLE 28|), modified-pole (BK 28l|). 



ISGW2 10|, and phase space (PHSP). Higher-resonance (and nonresonant) channels are used to 
match the sum of exclusive semileptonic branching fractions to the inclusive semileptonic branching 
fraction. 



Channel 



B (%) 



Form factor 



Comment 



DO 



K* 



2.16(17) [1] 
3.50(5) [5] 
0.11(11) 
0.11(11) 
0.12(3) [17, 29] 
0.288(9) [5] 
0.16(2) [2] 



SPOLE 
BK 

ISGW2 
ISGW2 
PHSP 
BK 

SPOLE 



rv = 1.62(8) and ra = 0.83(5) [17] 

aBK = 0.30(3) [5] 

B from Ref. [10] scaled by Ref. [5j 

B set to same as — > K^e'^Ve 

Nonresonant 

aBK = 0.21(7) [5j 

rv = 1.4(3) and rs = 0.6(2) [2] 



10 form-factor 



Higher-resonance decay branching fractions are predicted by the ISGW2 
model and remaining gaps are filled by nonresonant decays. We assume the size of the non- 
resonant component of D — )■ K7re~^i>e to be about 5% [o, 2^. Uncertainties of unobserved 
higher-resonance channels are assumed to be ±100% of the predicted branching fractions. 

The expected leptonic decay contributions due to the r+ — t- e~^h'ei^T decay in and Z)+ 
are used to correct nonsemileptonic electrons in our measurements as shown in Tables IVIIII 
and[IXl For decays, this component is expected to be large, and we include the leptonic 
decay component in the extrapolation fit. 
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TABLE VIII: Summary of semileptonic decays used to perform the spectrum extrapolation. 
Assumed branching fractions are shown in the second column; normalization of each component 
is allowed to float within the given uncertainty. Form-factor (FF) models used to describe the 
shape of each spectrum are shown in the third column: single-pole (SPOLE [28]), modified-pole 



(BK [28]), ISGW2 [1(|, and phase space (PHSP). Higher-resonance (and nonresonant) channels are 
used to match the sum of exclusive semileptonic branching fractions to the inclusive semileptonic 
branching fraction. The size of the expected secondary electron contribution from the leptonic 
decay — >■ t^i't is shown in the last row based on the known branching fraction of — )• /i"^!/^ 
decay 27] scaled by the standard model decay rate ratio [l3] T{D~^ — )• T^i'r)/T{D^ 



2.67. 


Channel 


B (%) 


Form factor 


Comment 


D+ - 


> if *Oe+z^e 


5.56(35) [3] 


SPOLE 


rv = 1.62(8) and rs = 0.83(5) [17] 


D+ - 


> K^e+i^e 


8.83(22) [5] 


BK 


asK = 0.30(2) [la 


D+ - 


> Kfe+i^e 


0.29(29) 


ISGW2 


B from Ref. [1^ scaled by Ref. [5] 


D+ - 


> Kfe+i^e 


0.29(29) 


ISGW2 


B set to same as Kfe'^Vg 


D+ - 


> Kne'^i'e 


0.32(8) [17, 29] 


PHSP 


Nonresonant 


D+ - 




0.405(18) [5] 


BK 


OBK = 0.21(7) [5] 


D+ - 




0.13(2) [4] 


BK 


FF set to same as — > TT^e~^iye [5] 


D+ - 




0.02(2) [4,i0, 30] 


BK 


FF set to same as — > TT^e~^iye [5] 


D+ - 




0.23(2) [2] 


SPOLE 


rv = 1.4(3) and = 0.6(2) [2] 


D+ - 




0.15(3) [2] 


SPOLE 


FF set to same as — > p^e^Uf, [2] 


D+ - 




e+i/ei^r 0.018 [17, 27] 







TABLE IX: Summary of Df leptonic and semileptonic decays used to perform the spectrum 
extrapolation. Assumed branching fractions are shown in the second column; normalization of each 
component is allowed to float within the given uncertainty during the fit. Form-factor models used 



to describe the shape of each spectrum are shown in the third column: single-pole (SPOLE [28|) 



and ISGW2 lOjj. The size of the expected secondary electron contribution from the leptonic decay 
Df T'^L'r is shown in the last row based on the known branching fraction of Df — > t'^i't 



decay 25|, |26|, and the shape is obtained from the evtgen jl6| ] MC program. 



Channel B (%) Form factor Comment 

D+ </)e+z/e 2.36(26) [T] SPOLE my = 2.1 GeV, tua = 2.28 GeV, 

rv = 1.849(112), and rs = 0.763(96) 
from Ref. [31] 

Df 7?e+i/e 2.48(32) [6] ISGW2 

Df r]'e+iye 0.91(33) [6] ISGW2 

Df K^e+Ue 0.37(10) [6] ISGW2 

Df K*^e+Ve 0.18(7) [6] ISGW2 

D+ ^ foe+Ue 0.40(6) [7, 32] SPOLE mpoie = 1-7 GeV [7] 

D+ T+i^r, T+ e+UeVr 1-003(79) [17, 25, 26] 
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